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Abstract 

Dilaton gravity with general brane localized interactions is investigated. Models with corrections 
up to arbitrary order in field derivatives are considered. Effective gravitational equations of motion 
at the brane are derived in the covariant approach. Dependence of such brane equations on the 
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of motion depends strongly on the symmetries assumed for the model and for the background. 
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1 Introduction 



Assuming that the space-time dimensionality is restricted to the number of currently known 4 
dimensions, the Einstein tensor appearing in the equation of motion of the standard theory of 
gravity is the most general rank 2, divergence-free symmetric tensor depending on the metric and 
its first and second derivatives only, while being linear in the latter. As for the presently testable, 
gravity-related phenomena, all experiments seem to show no discrepancies between the General 
Relativity and observations^]. Nevertheless, a beyond the Standard Model theory is commonly 
pronounced necessary, with string theories amongst the most serious candidates. As compared 
to the standard Einstein theory of gravity, the effective Lagrangians obtained within the string 
theories framework contain corrections with higher powers of the Riemann tensor [2]. Moreover, 
the a' expansion in string theories predicts corrections of the higher order in derivatives not only for 
the metric tensor, but also for other fields, such as the dilaton. In this work we investigate models 
with higher order corrections for both the metric and the dilaton, as well as mixed gravity-dilaton 
interactions. 

The order of the corrections is restricted by the dimensionality of the space-time. Specifically, 
the 7V-th power of the Riemann tensor can be included into a dilaton gravity Lagrangian only if 
the number of space-time dimensions, d, is sufficient, i.e. for d > 2N. For example, the 2nd order 
contribution quadratic in the Riemann tensor, known as the Gauss-Bonnet (GB) term [3J, can be 
taken into account in a dilaton gravity theory already in 4 dimensions. On the other hand, in 
the classical, no-dilaton theories, the GB term is a full divergence in 4 dimensions. It becomes 
dynamically relevant in an at least 5-dimensional space-time only. For pure gravity the GB term 
was generalized to higher orders^ by Lovelock [Bj. In a previous paper by the present authors [7j, the 
Einstein-Lovelock theory of gravity was generalized by coupling it to the dilaton. The appropriate 
action and equations of motion were constructed up to arbitrary order in derivatives of both the 
metric tensor and the dilaton. 

The main motivation for the dilaton gravity theory constructed in [7] comes from string theories, 
which are considered to be the most promising attempt to quantize gravity and unify it with all 
other known interactions. The gravitational sectors of the effective field theories derived from 
string theories coincide with the Einstein theory of gravity only in the lowest order. Beyond such 
approximation they involve higher powers of the Riemann tensor as well as interactions with other 
fields with more than two derivatives. The dilaton gravity model considered in the present paper 



1 However, there are claims that some alternative models could provide better explanation for the late-time 
cosmic acceleration and/or for the phenomena attributed usually to the existence of the Cold Dark Matter, (see e.g. 
[I] and references therein). 

2 Explicit formulae for the 3rd and 4th orders expressions are given in [3] and [5], respectively. 
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is closely connected to string theories - its interactions with up to four derivatives are exactly as 
those present in the effective theory derived from string theories and restricted to the gravity and 
the dilaton field (see e.g. [5]). Such correspondence has not been proven for interactions with six 
or more derivatives. However, there are considerable indications (see the discussion of the 0(d,d) 
symmetry in [JJ) that our model may be a part of the effective string dilaton gravity action also at 
the level of more than four derivatives. 

Dating back to the M-theory based work by Hofava and Witten [9], the idea of localizing 
the Standard Model on a brane embedded in a higher- dimensional space-time [10] gained quite a 
lot of attention. Corrections to the gravity interactions at the brane due to the bulk fields were 
investigated by many authors [11]. It seems interesting and important to consider what gravity 
would be induced at the brane, if the bulk action was given by a higher order dilaton gravity theory. 
Hence, the purpose of this work is to derive effective gravitational equations at the brane for the 
models of arbitrary higher order in derivatives, constructed in [TJ. In order to keep full generality, 
the procedure will be carried out in the covariant approach. Starting from the rf-dimensional (d > 4) 
higher order dilaton gravity theory, equations of motion in the effective (d — 1) -dimensional theory, 
i.e. for a co-dimension 1 brane, will be derived. 

For the standard lowest order gravity the effective equations of motion at the brane were derived 
in the covariant approach in [12J. That analysis was extended to the GB gravity in [13]. The 
effective equations at the brane for the lowest order dilaton gravity were derived in [T3] (however, 
not in a fully covariant way). The covariant approach was employed in JT5] for cosmological 
applications of dilaton gravity. Certain second order gravity models with branes, including either 
first or second order corrections for the scalar field, were analyzed for specific metrics in e.g. [16] 
and [17], respectively. Although brane models for arbitrary order Einstein-Lovelock gravity were 
investigated in [T5], neither the dilaton field was included, nor the covariant approach was adopted. 
No comprehensive results for theories which simultaneously take into account interactions of the 
higher order in the Riemann tensor and involve the dilaton have been presented so far. 

The main goal of the present work is to obtain the effective brane equations of motion for 
arbitrary order dilaton gravity models using the covariant approach. Due to the already pointed 
out relation of the higher order dilaton gravity model [7J to string theories, the covariant derivation 
of the effective brane equations in this setup will allow for studying in a more systematic way the 
potentially observable effects of string theories e.g. on specific cosmological models. 

The rest of this paper is organized as follows: In section [2] we define the d- dimensional dilaton 
gravity models with corrections of the higher order in derivatives and with brane interactions whose 
exact nature we choose not to specify. The bulk equations of motion for these models are derived 
and rewritten in terms of quantities either projected on the brane, or perpendicular to it. Section 



3 



[3] is devoted to the analysis of the junction conditions at the brane. The effective equations of 
motion at the brane for a general case with arbitrary order corrections are obtained in section 
HI Construction of all possible brane equations is carefully discussed. The importance of the 
bulk Z2 symmetry for the form and existence of the effective equations is analyzed. In section [5] 
explicit results for the lowest order dilaton gravity and for the theory with up to four derivatives 
are presented. Some of the results for the latter, because of their complexity, are moved to the 
appendix. Section [6] contains our conclusions. 



2 Higher dimensional equations of motion 



We consider the d- dimensional dilaton gravity theory described by the following Lagrangian: 

N< max 
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This is the Lagrangian (in the notation explained below) of the higher order dilaton gravity con- 
structed in [7] and generalized by including a bulk scalar potential V(4>) and general brane inter- 
actions given by Lb- The position of the branej is described by the Dirac delta type distribution 
5b- As was discussed in the Introduction, the main motivation for considering this dilaton gravity 
model is its close relation to string theories. Investigation of its properties in a brane scenario is a 
step towards studying e.g. the potentially observable cosmological effects of string theories. 

The above Lagrangian is written in a very compact form using the following notation introduced 
in [TJ: T is a generalization of the ordinary trace. Acting on an arbitrary rank (m, m) tensor M it 
returns a number given by 

t(m) = 5;i%:::%M^--^ ai(T2 ... am , (2) 

where 5 with m pairs of indices is the generalized Kronecker delta 
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Asterisks are used as indices in the Lagrangian (JT]) to indicate ranks of tensors (e.g. to distinguish 
the Riemann tensor from the Ricci tensor or the curvature scalar). Under the generalized trace T 
there are powers of "sums" (denoted by ©) - i.e. of linear combinations of different ranks tensors. 



3 The theory can include more branes, but we choose to focus our considerations on the brane for which we will 
be deriving the effective gravitational equations of motion. 
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This should be understood as a compact notation for the following operation: First, we perform 
the algebraic manipulations (sums and powers), treating all tensors as ordinary numbers. Then the 
generalized trace ([2]) of the obtained linear combination of tensors products should be understood 
as the corresponding linear combination of the generalized traces of tensors products. For example: 



T 



1 -Tp ** 

2 '^-** 



2(w): 



ir(7e:;?e;:) + 2r(7e::(vv)» +4T((vv):0(vv): 



1 rciCT20-3CT4 <1->piP2 

4 P1P2P3P4 <71<T2 /V 0"30"4 



I" ^ U P1P2P3 /V O"!^ 



where, due to the properties of the Riemann tensor and the second covariant derivative of a scalar, 
we can unambiguously define: 1l p ° = Vf a ^ and (VV)^0 = V v d^4>- 

With the above formulae it is easy to check that the number of higher order terms in the 
Lagrangian ([1]) is not arbitrary. The term in <^f with the generalized trace ([2]) of the highest rank 
tensor is proportional to 
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It is obvious from the definition ([3]) that, because of the antisymmetry in all indices of one type 
(covariant and contravariant), the (2iV max , 2iV max ) rank generalized Kronecker delta in eq. (jl]) 
is non-zero only if the number of space-time dimensions d is sufficiently large. Thus, there is an 
upper limit on the order of the interaction terms present in our Lagrangian: iV max < (d/2). 

Using the results of [TJ, bulk equations of motion can be derived from the Lagrangian (CQ). They 
read 



■Afi max 



N=l 

N„ 



i7C©2(VV):0©(-l)(<90) 5 



N 







N=l 



A? 



tJ>b = 



(5) 



(6) 



where the brane localized terms, Tn U and ta, are calculated from the brane Lagrangian namely 
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In the tensor equation of motion (JSJ) we introduced another generalization of the standard trace 
given by 

%{M) = S^^M^^ nn .^ m . (8) 
It maps arbitrary rank (m,m) tensors into rank (1, 1) onea_. 



For m = d the generalized Kronecker delta fiuplpZ'-'-'-p™ should be replaced in ([5]) with the following combination: 
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The main goal of the present work is to find the effective {d — l)-dimensional equations of 
motion at the brane located at the support of 8b- The first step in this direction is to identify 
parts "parallel" and "perpendicular" to the brane in all relevant tensors. In order to achieve this 
in a covariant way, we start with introducing a vector field n M normalized to 1 and perpendicular 
to the brane at its position. The choice of is not unique, due to the freedom in the bulk and 
the sign ambiguity at the brane. However, as we will see later, the effective brane equations of 
motion are unique. With any such a vector field we define the metric = g^ v — n^n v . This 
expression holds throughout the d- dimensional space-time but, restricted to the brane position, it 
just yields the metric induced on the brane. Subsequently, we divide the o?-dimensional tensors into 
parts parallel and perpendicular to the vector field n p as 

K°p = R 7p~ 2 K 1» K p] - 4 nlVD l»K] ~ 4 n l» DiU K] ~ 4 n[l/ n { ,£ n K;l + 4 nK [lM (KK);\ , (9) 
V" V M = D v D lt <f> + K;£ n cj> + n v D^£ n <P + n M D"£ n <f> - n u K°D a (j) - K v a D°$ 

+ nX(A 2 0-a A V A 0), (10) 

m 2 = w) 2 + (^) 2 , (ii) 

where K^. = \£ n h iiv is the extrinsic curvature of hypersurfaces orthogonal to n M , £ n is the Lie 
derivatives along n M , and a x = n p V p n A . It is important to distinguish between the d- dimensional 
tensors associated with the full metric g^, namely the Riemann tensor TZ^ and the covariant 
derivative V M , and the (d — l)-dimensional tensors associated with the metric h^ v restricted to the 
brane position, i.e. and D^. 

The right hand sides of eqs. (l9l-(jTTT) are expressed almost entirely in terms of the vector field 
ra^, Lie derivatives along and brane quantities orthogonal to n/ 1 , i.e. Rp£, and D^. The term 
a x V\4> = n p (V p n x ) (V \<p) , containing the d- dimensional covariant derivatives, is the only exception. 
However, as we will see later, the derivation procedure for the effective gravitational equations at 
the brane will be constructed in such a way, that a x V\(f) will not appear in our final results. 

The Riemann tensor IZ 1 ^ and the tensor of the second covariant derivative of the dilaton V u V fJi <f) 
appear in the bulk equations of motion ([5]) and ((6]) exclusively under the generalized traces ([2]) and 
OH]), which involve full anti-symmetrization in all covariant and contravariant indices. Hence, the 
projection equations OH]) and ffHJ]) can be simplified when used under those traces, namely 

k: -+ R::-2K:K:-4(nn):(£ n K:-(KK):)-8(nD):K:, (12) 

(VV):0 [(DD):<p + K:£ n <p}+(nn):(£ n 2 <p-a x V x <p)+2[(nD):£ n <p-(nKD):<p], (13) 



5 For a given tensor M^lal-.-aT an( i arbitrary direction u M , the Lie derivative along is denned as follows: 

£ v M£g:::sr = v x v x m£%::z? - YZx m&£:±::StVav« + xj =1 ai^z^T v ct v x , thus e.g. £ n <t> = n A v A ^. 
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where in order to make the formulae more compact we introduced the following notation: (nn)* = 
n*n*, (VV): = V*V*, (DD)l = D*D*, (KK)t = K* X K$, (nD)l = \ (n*D* + n*D*) and (nKD)l = 
\{nJ<lD x + n*K^D x ). 

Employing the dimensional reduction formulae (JT2|) . (|T3|) and (fTTl) . the bulk equations of motion 
(El) and (El) can be rewritten as 
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These equations remain defined and valid throughout the full dimensional space-time. They are 
smooth in the bulk, but contain discontinuous and singular (distribution-like) terms localized at 
the brane. The effective brane equations of motion can be obtained if the information from both 
smooth and singular parts of (fl4|) and (1151) is properly taken into account at the position of the 
brane. 



3 Junction conditions 

In order to find the junction conditions which have to be fulfilled at the brane, we shall define 
a 1-dimensional integration in the direction perpendicular to the brane. In order to achieve this, 
we will employ a special family of integral curves associated with the vector field n M . Specifically, 
let us define the curves 7 M (A) by requiring the vector field to be at each point tangent to 7 M : 
(d7 At /dA)(A) = n M (7 i/ (A)). We choose the parameterization of these curves in such a way that each 
of them crosses the brane at the point 7 M (0). For a scalar field f{x^), at each point we consider 
the following integral: 

/ /(*")= / A2 /(7£(A))dA, (16) 
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where T^(Ai, A2) = {7z(A)|A G [Ai, A2]} is a part of the curve 7^ (A) crossing the brane at the point 



x M . This integration is "inverse" to the Lie derivative^ in the following sense: 

A/K)=/(7^(A 2 ))-/( 7 ^(A 1 )). (17) 

This integral can be generalized to the case of arbitrary tensor integrands by applying appropriate 
pullback and pushforward operations along the family of curves 7 M (A). 

Junction conditions for a given point Xq on the brane are obtained by integrating the d- 
dimensional equations of motion along (an infinitesimal part of) the curve 7^ (A) going through 



•'V 



lim 

e-S-0 



(18) 

J-B " u Jr% (-e,+e) 

where F(x M ) is the scalar or the tensor given by the left hand sides of the bulk equations of motion 
flUD and (USD, respectiveljfl. 

Only certain terms in the (i-dimensional equations of motions will not yield zero during the 
above infinitesimal across-the-brane integration, namely the brane contributions involving r^ v and 
r<f,, explicitly proportional to 5b, and the terms containing Lie derivatives (i.e. £ n 2 (J) and £ n K lJbV ) 
of quantities which can be discontinuous at the brane (£ n <f) and K^, respectively). Hence, the 
relevant parts of the bulk equations of motions (I14p and (I15p yield 



N- max f* 

N=l J ±B 



nn 



)l£n{K 



■1)£, 



1 D** 

2 ** 



Nmax p / 

y2a N N r[£ n (Kt 

N=l J ±B V 



2 ** 



2( J DD):0 © (-1) W) 2 j © (-1){K © {-1)£, 

~l)£n<£) 

2( J DD):0 © (-1) (D<f)) 2 ) © (-i)(/C a 



N-l 



T^u , (19) 



-D£r 



N-l 



T+. (20) 



The integration of the terms containing the brane localized sources t^ u and was easy due to the 
obvious (defining) properties of the Dirac delta distribution 5b- The remaining integrations require 
a more careful treatment. The integrands involve products of the distribution-like objects, £ n K flv 
or £ n 2 (j), and the potentially discontinuous objects, K^ v or £ n (j). Strictly speaking, such integrals 



6 From now on whenever we refer to a Lie derivative, it should be understood as the Lie derivative along the 
vector field n M . 

7 Integral (|18|) is well defined for any point of the space-time. However, away from the brane the integrands arc 
smooth and the integral vanishes in the e —> limit. Non-trivial results are obtained only on the brane, where there 
are localized sources and r^. 



have no mathematically unambiguous meaning, as distributions are defined by their integrals with 
smooth functions. Therefore, we need some kind of regularization to deal with such terms. Fortu- 
nately, there is an obvious way to regularize and calculate the integrals in eqs. ( II 9p and ( I2(jp . Let 
us first consider terms of the form sf k £ n f, where s is smooth, whereas / may be discontinuous at 
the brane. Employing the formula (1T7|) and the Leibniz rule for the Lie derivative we get 

/ sf£ n f = ^ f s£ n =&[ £ n (sf k+1 ) = j^[sf k+1 U • (21) 

J±B J±B J ±B 

The second equality follows from the fact that the integral of (£ n s)f k+l vanishes in the e — > 
limit. The square bracket with the ± subscript was introduced to denote a jump in the value of 
some quantity when crossing the brane: 

UK)}± = [/«)]+ - lf(<)}- , (22) 

while the square brackets with the subscripts + and — denote the limits of a given bulk quantity 
when approaching the brane from the "+" and the "— " sides, respectively: 

[/«)]+(-)= lim /( 7 £ (e)). (23) 

We intend to apply the regularization (12 ip to the integrals in eqs. (|T9|) and (120]) . It is instructive 
to start from considering a simple example of such a calculation: 
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The same result is obtained if the regularization (I2~TI) is used when treating the formal sums of 
tensors of different ranks, under the generalized traces T and l~^ u , as ordinary functions. Performing 
similar calculations for all terms present in ( JT9l) and ( l20l) . we obtain the junction conditions as 
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These formulae can be rewritten in a slightly (but qualitatively) different form in the case of theories 
with the bulk Z 2 symmetry and the brane located at the Z 2 symmetry fixed point. Then, for any 
Z 2 -odd quantity /, [/]_ = — [/]+ and [f}± can be replaced with 2 [/] + . 

A comment on distinguishing both sides of the brane is in order here. We started our construc- 
tion with the vector field n M , which at the brane is normal to it. The vector field — n M has the same 
property. Replacing n M with — n M (which corresponds to interchanging the "+" and the "— " sides 
of the brane) we have to change the sign of A in the family of curves 7^ (A) used in (1TB]) to define 
the integral j ±B . From (1221) it follows that this in turn changes the sign of [■ ■ - ]±. The vector 
field n M enters linearly also the definitions of the Lie derivative £ n and the extrinsic curvature K^. 
Hence, the combination (K* © (—1) <£„</>) changes sign together with n M . Observe that only odd 
powers of this combination appear in the junction conditions (j24p and (1251) . Thus, the change of 
sign in the definition ( I2"2"j) of [■ • • ]± is compensated by the change of sign of £ n <j) and K^ v . The left 
hand sides of (j24j) and (I25p do not depend on the sign of n 11 - as it should be, as the brane localized 
interactions (JTj) appearing on the right hand sides of these equations do not depend on which side 
of the brane is called the "+" one. 

The junction conditions (1241) and (1251) determine the jumps in the values of the extrinsic curva- 
ture, K^v, and the Lie derivative of the scalar field, £ n <f), caused by the brane localized interactions 
described by and r^. Unfortunately, the solutions for [if^i and [£ n <f)]± (or [i^ AM/ ]+ and [£ n 4>] + 
in the case of Z 2 -symmetric models) can be found explicitly only in some simple cases, for example 
when iV max = 1 or when we are considering solutions which are highly symmetric. 



4 Effective brane equations of motion 
4.1 Definitions 

In the previous section we found the junction conditions (124"]) and (1251 which have to be fulfilled at 
the brane. Presently we would like to obtain (d — 1) -dimensional equations which we could call the 
"effective brane equations of motion" . First of all, a precise definition of such effective equations 
of motion at the brane has to be given. There are two obvious properties of these equations: First, 
they must follow from the full d-dimensional equations of motion ( TT4]) and ( TL5"]) . Second, they should 
describe the behavior of the quantities defined exactly on the brane, or infinitesimally close to it. 
We can obtain such equations in a very similar way to that employed in the junction conditions 
derivation. Specifically, for each point Xq on the brane, we integrate the full d- dimensional equations 
of motion (JT3j) and (TToT) over some infinitesimal interval "perpendicular" to the brane. However, the 
integration (jT6l) has to be generalized in order to obtain something different than just the junction 
conditions (124]) and (125]) . An obvious way to generalize the integral ([16]) is to use a weight function 
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w(X): 

V2 f(^(X))w(X)d\, (26) 



Ai 



where 7^ (A) is again the integral curve of the vector field n M intersecting the brane at the point 
Xq. In fact, we need a family of weight functions w e (X) such that the support of w £ (X) is included 
in the interval (— e, +e). Any such a smooth function can be written as the following sum: 

We (\) = w p + wi~\X) + wi + \X) , (27) 

where w| ' is a constant (finite in the e — > limit), w[ \X) has the support in (— e, 0) and wi + \X) 
has the support in (0, +e). The brane equations of motion are then obtained by integrating the 
d- dimensional equations of motion (|T4|) . (1151) like in (|26|) and taking the limit of e — > 0, namely 

lim / F( 7 £ (A) W(A) dA = Co / F(a#)+c_ lim F(x 1 +er^)+c+ lim F(a#+en"), (28) 

where the coefficients c depend on the chosen w e : 

^ _) (A)dA, c+=lim/ w e (+) (A)dA. (29) 

- e Jo 

By adopting different w e we obtain three independent (d — 1) -dimensional equations. The first 
term on the r.h.s. of ( 128|) is proportional to the corresponding junction condition discussed already 
in the previous section. The remaining two terms are proportional to two directional limits of the 
bulk equations (i.e. the limits of the bulk equations when approaching the brane from the "+" or 
the "-" side). 

None of these equations can be called an effective brane equation of motion. The directional 
limits of the bulk equations have no explicit dependence on the brane localized quantities ([7]). The 
junction conditions do depend on the brane localized quantities, like the brane energy momentum 
tensor t^ u , but they determine the jumps in the bulk quantities values and not the dynamics of 
the brane quantities. However, one should not draw the conclusion that it is not possible to define 
any effective brane equations of motion, as there are more <i-dimensional bulk equations than are 
needed for a (d — l)-dimensional brane gravity. In addition, one should remember that the metric 
tensor and the scalar field are continuous at the brane. The effective equations of motion at the 
brane are obtained by combining all available information: the junction conditions, the directional 
limits of all bulk equations and the continuity conditions. This will be explained in more detail in 
the next subsection - employing our model with the bulk equations of motion (TH|) and f|T5|) . 

4.2 Construction — general case 

We start with projecting the bulk tensor equation of motion ( Tl4|) on the brane hypersurface and/or 
on the normal vector field . Multiplying f|T^|) by h^h^, using the fact that n M is orthogonal to all 
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brane directions, and taking appropriate limits of the bulk fields, we get the following two equations 
(one for each side of the brane): 
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where we introduced symbols A4 and M to denote the following formal combinations of the tensors 
of different ranks: 



M = \RZ © 2(DD):<f) © (-l)W) 2 © (-1) {K © (-l)£ n <PY , 
AT = [D*K: © (K$D X (f> - D t £ n 4>)] [D*K © (K* x D x </> - D*£ n 



(31) 
(32) 



which will appear under the generalized traces T and in numerous formulae. The index 
structure of M is different than that of all other objects appearing under T and T^ v . So far all 
such objects had natural pairs (covariant-contravariant) of indices - corresponding to the pairs of 
indices present in the generalized Kronecker delta ([3]). In TV there is one unpaired covariant index in 
the first square bracket in (1321) and one unpaired contravariant index in the second square bracket. 
It should be understood that these two indices form a pair under T and T^ v . Thus, we have for 
example 

TW = &;iZl i D ^ K ? 2 ) {D P1 KD + iP^K") (K p x 1 D x <p - D pi £ n (f)) +... 

The bulk dilaton equation of motion ( Tl"5]) is very similar to the projection ( 130]) of the bulk tensor 
equation of motion on the brane (i.e. obtained by contracting with h^h v a ) and reads 



N=l 



2N 



{(KK): - £ n Kl) © (£ n 2 <j) - a x V x 
N(l -N)MM N ~ 2 



M 



N-l 



0. 



(33) 



-(+) 



The contractions of the bulk tensor equation of motion (Tl4l) with h^n u and n^n u are given by, 
respectively, 



N=l 



£ a N T pa ( N D°K: © (KZD x <j) - D a £ n 

N max 



N 



N=l 



-(+) 



-(+) 







0. 



(34) 
(35) 
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Before discussing the properties of the above equations of motion, let us rewrite the tensor 
equation fl30l) in a somewhat different form. One of the reasons is that we would like to remove 
from our equations all terms containing the quantity o a Va0, which is the only term being neither 
"parallel", nor "perpendicular" to the brane. Another reason is that we would like to rewrite the 
tensor equation fl30|) in a form more suitable to compare our results with those presented so far 
in the literature. In order to achieve these goals, we conduct the following procedure: Using the 
decomposition (jHJ) and the definition of the <i-dimensional Weyl tensor C^ upa : 

Cpupa = Ttfiupa — {9 n\f>R'<T\v ~ 9v\p^a\p^ + (d-\)(d-2) 9p[p9a]vTi i (36) 

we obtain the equation 

£ n K iiV - (KK)^ = ^ V [(h£ n K) - (KK)] - ^ [R^ - KK^ + (KK)^} 

+ V [R~K 2 + (KK)} - ^E, v , (37) 

expressing the Lie derivative of the extrinsic curvature, £ n K fJiU , as a function of its trace, (h£ n K) = 
h tJ,v £ n K fJiV , and the following projection of the d- dimensional Weyl tensor: 

= n a h^hl C aM& . (38) 

After using eq. (]3~T]) , the tensor equation fl30|) depends (explicitly) only on two second Lie derivatives 
of the bulk fields, (h£ n K) and £ n 2 (j), which appear always in the combinations [(h£ n K) — (KK)] 
and [£ n 2 (j) — a x \7 \(fi]. These two combinations can be calculated from a system of two linear 
algebraic equations consisting of the scalar equation (13"3"|) and the trace of the tensor equation ( 13"Uj) 
- both evaluated "next to the brane" . The solution of this system of equations reads 

(h£ n K)-(KK)=(d-l) b f- b lf\ (39) 

f> 2 i An J, biBi—b 2 B 

£ n - a V A = — T2— , (40) 

b b 2 - b{ 



where 



Nm ax 

b m = Y,a N NT((h:) m M N - 1 ), (41) 



jV=l 
N m ax 



JV=1 



Bm = J2^fr((K) m M N ®^PM N - 1 ®8N(l-N)AfM 



N-2 



-(d-irV(<f>) + (l-m)V'(<f>), (42) 
V=K~ KK + (KK)l - [R-K 2 + (KK)] + (d - 2)K , (43) 
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whereas M and M are defined in (13T|) and (|32|) . respectively. Substituting eqs. (1371) . ( )39|) and ( 140]) 
into the projected tensor equation f )30|) . we obtain a tensor equation of motion without explicit 
dependence on the second Lie derivatives of the bulk fields, namely 



N max 

« pa v 2 9 1 

N=l 



fj i/ 



M 



N 



N-l 



2N 



h 



«6 igo-6ogi 
b b 2 - b\ 



IN(1 - N)MM N ~ 2 
b x B x - b 2 B 



b b 2 - b\ 



M 



N-l 



(44) 



-(+) 



It is useful to reformulate this equation further in order to rewrite it as a sum of the ordinary 
lowest order Einstein equation and some corrections. The first step is very simple, we just have 
to multiply (jUJ) by (d — 3)/(d — 2) in order to get the usual normalization of the Ricci tensor. A 
second step is necessary, as the ratio of the coefficients in front of h pa R and R pa is different than 
the desired —1/2. To solve this problem we have to add a product of eq. (13"5"|) and the brane metric 
tensor h pu with the appropriate coefficient. Finally, we obtain the following Einstein-like tensor 
equation of motion: 



(d-3)(2d-3) , y 

(d-l)(<f-2) ,b P° V 



d=s — /i^r 

N=l 



2N 



MB 



E 

Af=l 

b B 1 



M 



N 



2N 
d-3 



VM 



b b 2 - b\ 



b x B x - b 2 B 
b b 2 - b\ 



A'n 



IN(1 - N) MM 



N-2 



N 



N=l 



0. (45) 



-(+) 



The above equation seems to be rather complicated, but, in fact, after calculating the generalized 
traces ([2]) and ([8]) it takes the most Einstein-like form we could obtain, as we shall show explicitly 
for iV max = 1,2. Moreover, it does not depend on the Lie derivative of the extrinsic curvature, 
£ n K^, and on the second Lie derivative of the dilaton, £ n 2 (j). The whole dependence on these 
quantities, which cannot be restricted by the junction conditions, is encoded in the projection ( 138]) 
of the bulk Weyl tensor on the brane, E pv . This dependence is implicit via the quantities B m and 
V defined in ( 14"2"|) and (|4"3"|) . respectively. Due to the explicit bulk contribution, represented by E^, 
the Einstein-like equation of motion (145 p does not form a closed system. Consequently, to describe 
fully the brane dynamics, the bulk solution usually would have to be known. 

All the equations we wrote down so far in this subsection are directional limits of the bulk 
equations of motion (or their combinations). In order to get any effective brane equations of 
motion we need some dependence on the brane sources ([7]). Such a dependence can be introduced 
by taking into account the junction conditions, which determine the jumps in the values of K pv and 
£ n (j), i.e. [-K" M i,]± and [<£ n 0]±, in terms of the sources T pu and r^. Let us now discuss what equations 
are obtained after substituting the solutions (not always known in a closed analytic form) of the 
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junction conditions flMD and ([25]) into the bulk equations ([30}, ([33]), ([34]) and (1331FI. 

The first obvious observation is that eqs. (130]) and ([33]) give no useful effective brane equations 
of motion, as they involve second derivatives of the brane metric and the dilaton field along the 
vector field n M normal to the brane, i.e. £ n K iiV and £ n 2 <fi, respectively. To show this more explicitly 
we consider the following construction: We choose some brane sources t^ v and r^, together with a 
metric h^ v and a scalar field <fi on the brane. More precisely, we choose [71^]+ and [0] + (which is the 
same as [/i^]- and [0]-, as the fields are continuous) at one of the two hypersurfaces infinitesimally 
close to the brane. Due to the junction conditions, the brane sources r^ v and restrict (or in some 
cases determine) the first derivatives in the direction perpendicular to the brane, i.e. + and 
[,£„</>]+ (or [-K'jui/]- an d [£ n (j)\-) at the chosen hypersurface "next to the brane". The values of these 
fields (i.e. [^^1+ and [0]+) and their first derivatives constitute simply the boundary conditions 
for the quasi-linear second order differential eqs. f [3"U]) and (13"3"]) . The Cauchy-Kowalewski theorem 
tells us that such boundary conditions problem can be solved at least in some neighborhood of 
the brane. The existence of the solutions does not require any relations between the brane fields 
and the brane sources. Without any additional assumptions (e.g. about the symmetries of the 
bulk/brane solutions) the eqs. ( 13U]) and (133]) and the junction conditions ( 12~4"]) and (123]) do not yield 
any effective brane equations of motion, as they are not sufficient to get any constraints on the 
brane fields and in terms of the brane sources t^ v and t$. 

From the above reasoning it is clear that in order to get any effective (d — l)-dimensional 
brane equations we need bulk equations of motion which do not involve the second derivatives 
perpendicular to the brane. Such equations, obtained from the d- dimensional tensor equation of 
motion by contracting at least one of its indices with the vector field n fl normal to the brane, are 
given by ( 134"]) and ([35]) . Let us start with the latter. The "+" side eq. ( [35]) is just a consistency 
condition on the "+" side first Lie derivatives of the brane fields, and [£ n (f>} + . Using the 

relations 

[Kfa,]- = [K^] + - [i^^]± , [£n<t>]- = [£n4>]+ - [£n<j>]± , 

the "— " side eq. ([35]) becomes another consistency condition on and [£ n <f)] + , this time 

depending on the brane sources t^ u and r^. It can be seen that in general, even in the simplest 
-^max = 1 case, these two consistency conditions obtained from eq. ( 133]) can be fulfilled, thus eq. 
( 133]) does not lead to any effective brane equations of motion. The situation may change if some 
assumptions about the bulk are made, e.g. if the bulk Z 2 symmetry, relating [^^1+ with [K^]^, 
and thus with and the brane sources t^ u and t$, is assumed. Such models will be discussed 

in the next subsection. 

We have shown that in general eqs. ( 130]) . ([33]) and ([33]) cannot be used to obtain any effective 
8 For this discussion it is better to use eqs. ([50]). (pt5|) and ([55]) . than the derived from them eqs. (pPT| and P51) . 
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brane equations. What remains to discuss now, are two eqs. (I34|) - one for each side of the brane. 
They were obtained from the full ci-dimensional tensor equation of motion ( fT4"j) by contracting one 
of the indices with and the other one with h^ u . The "+" side eq. (I34p gives (d — 1) additional 
consistency conditions on [i^" A1I/ ] + and [£ n (f)] + . Together with those 2 resulting from eq. (|35|) . we 
have (d + 1) such conditions. Hence, the number of the conditions is smaller than the number of 
the (gauge independent) degrees of freedom in and [£ n (p} + , i.e. these conditions can be in 

general fulfilled. 

The situation is more interesting for the "— " side eq. ( I34"l) . One could try to use the same 
argument as for the "— " side eq. ( 135|) and expect to get additional consistency conditions for + 
and [£ n <f)] + . However, this time that argument does not work. The reason is that the difference 
between the "+" and the "— " sides projections f)34p is closely related to the junction conditions 
( 12"4"|) and fl2^|) . After a somewhat tedious calculation, it can be shown that such a difference is 
equivalent to the following condition on the brane sources: 

D^ u + (V"^ - vO = • ( 46 ) 

This condition has the same simple form for all theories of the structure described by the Lagrangian 
([1]) of arbitrary higher order in derivatives. It does relate the brane sources ([7]) to the brane metric 
(the covariant derivative on the l.h.s. is covariant with respect to the induced brane metric h^ u ). 
However, its character depends crucially on the brane localized contribution £g to the Lagrangian 
([1]). For a wide class of Lb eq. (146]) is not a dynamical brane equation of motion, because it 
does not involve the second derivatives (in the brane directions) of the brane fields. It is rather a 



consistency condition on the source termed t^ u and r^. Such consistency conditions are typical of 
gravity theories |19j . 

The situation changes for more complicated brane Lagrangians Lb- Specifically, if Lb contains 
e.g. localized kinetic terms for the gravity [20] or for the scalar [21], the condition ( H6|) becomes a 
dynamical brane equation of motion involving the second derivatives (in the brane directions) of 
the brane fields h^ u and <fi. Certainly even in such a case we do not get a full system of the brane 
gravitational equations of motion. The reason is obvious: the condition (H6|) is a (d — l)-dimensional 
vector of equations, while the brane gravitational equations should have a (d— l)-dimensional tensor 
character. 

The condition (T4"6"j) seems to be of the lowest order, as all the higher order terms present in the 
bulk Lagrangian (Cp) "canceled-out" in its derivation when we combined the higher order junction 
conditions ( 124)) and (|25|) with the higher order bulk equation ( l34j) . However, the condition (146]) may 



9 In the case without the dilaton, i.e. in the case of the Einstein-Lovelock theory of gravity, eq. reduces to 
the requirement that the brane tensor source t^ v must be covariantly conserved. 
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involve terms with more than two derivatives of the fields if appropriate interactions are present in 
the brane Lagrangian C B . 

Let us summarize the general case without Z 2 or any other symmetry imposed on the bulk 
Lagrangian ([1]) or the bulk background solution. We have the following system of equations: 

(i) Junction conditions and (125j) . which determine (not always explicitly) the jumps in the 
values of the extrinsic curvature and the Lie derivative of the dilaton: [if^i and [£ n 4>]±- 

(ii) System of equations for and [£ n 0] + given by the "+" side eq. ( IMj) together with both 
eqs. (135p . The number of these equations is smaller than the number of (gauge-independent) 
degrees of freedom in [1^]+ and [£ n <J)] + , so that system of equations in general has solutions. 

(iii) Consistency condition (116]) on the brane sources r Mi , and r^. 

The first two systems of equations, (i) and (ii), are not restrictive enough to obtain any effective 
brane equations of motion relating the dynamics of the brane fields to the brane sources (j7|). One 
can choose any brane fields, i.e. induced metric and dilaton </>, and any brane sources, and T<j>, 
satisfying the consistency condition (iii). In general, there exist values of [K fJiU }_^ and [£ n 0]_( + ) 
which fulfill eqs. (i) and (ii) for any such choice. On the other hand, the consistency condition (iii) 
may have a character of a dynamical brane equation of motion. This depends crucially on the form 
of the brane localized interactions described by 

4.3 Construction with Z2 symmetry 

The bulk Z 2 symmetry is employed in many papers on brane models. We will show now that such 
a symmetry not only simplifies the calculations, but can also change qualitatively the problem of 
the existence of the effective brane equations of motion. 

The bulk Z 2 symmetry relates the fields at the "+" and the "— " sides of the brane. We have 

e.g. 

[K^U = = \ [K^} ± , [£ n <f>] + = -[£ n <f>]- = \ [£ n <P}± ■ (47) 

In such a case, the "— " side eqs. fl30l) . (l33]) - (]35]) coincide with the "+" side ones. Moreover, the 
junction conditions ([21]) and (l25j) determine the extrinsic curvature K^ v and the Lie derivative of 
the dilaton £ n <f) on both sides of the brane. The junction conditions are given by equations of 
order (2iV max — 1), so in general they cannot be solved explicitly for models of the higher order 
in derivatives. They can be solved analytically in some cases of highly symmetric configurations. 
Otherwise we have to solve them numerically. Substituting the (explicit or not) solutions of the 
junction conditions, [K flu ] + (T fJiU , r^) and [£ n (f)] + (r^ u , r^), into eqs. ( |35]) we obtain the following 
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effective brane equation of motion: 



N=l 



±ic© 2{DD)i<j> ©(-i)w) 2 



1) [/C] + ( V ,r )©(-l)[£ n ^ 







(48) 



This is the first equation we found which has the character of an effective brane equation of motion 
even for simple brane Lagrangians Lb- It involves the brane sources (J7|), as well as the second 
derivatives (along the brane directions) of the brane metric, an d the dilaton, i.e. D v D ll (\) and 
(D<ft) 2 . The covariant derivatives of the dilaton can be eliminated using the condition (|46p . This 
way we obtain one equation relating the dynamics of the brane metric h^ u to the brane sources t^ u 
and r,f,. It is the only such equation which appears in our theory - if we assume nothing else but 
the Z 2 symmetry in the bulk. 

One effective equation of motion in a (d — l)-dimensional brane gravity is not much. The 
induced brane metric h^ v has d(d — l)/2 gauge independent degrees of freedom. Nevertheless, 
even one equation of motion can be very important if we restrict our attention to highly symmetric 
brane solutions. For example, any maximally symmetric space-time (like de Sitter space-time used 
to describe inflation) is fully determined by just one parameter - the curvature scalar. Similarly, 
the cosmologically important Friedmann- Robert son- Walker space-time depends on one function of 
time only - the cosmic scale factor, for each sign of the brane spatial curvature. 

Despite the fact that there is only one "true" effective brane equation ( 148]) . it is sometimes 
convenient to write down the full tensor Einstein-like equation of motion for the brane fields h^ v 
and 0. It follows straightforwardly from eq. ( 145]) and is of the same structure, namely 



(rf-3)(2d-3) 



h oa V 



(d-l)(d-2) '*P<r 



d-2 2 V 

N=l 



M 



N 



2N 



h 



b 1 B -b B 1 bJh-b^Bp 
b b 2 - b\ b b 2 - b\ 



M 



N-l 



d- 3 
d- 



8N(1 - N)MM N ~ 2 
a N 



fE 

N=l 



h pa Tl. 



N 







(49) 



with the metric h^ v and the dilaton <fi taken on the brane and all (implicit) terms involving [A' Miy ] + 
and [£ n 0] + (with the bulk Z 2 symmetry it does not matter which side of the brane we choose) 
replaced with the appropriate solutions of the junction conditions (1241) and (125]) . Obviously, the 
above effective brane equation does not constitute a closed system. The solution of the equations 
of motion for the bulk gravity has to be known to fully describe the gravity induced on the brane. 
However, the total bulk dependence of the above equation is described by the projected Weyl tensor 
E^v, given by ( 138]) . which enters (implicitly - via V and B m defined in ( 143]) and ( 142]) . respectively) 
in a quite complicated and in general non-linear way. 
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Our effective brane equation ([49]) has one feature which is unusual for the standard equations 
of motion, but quite typical of the brane ones. Specifically, there are no terms linear in the brane 
energy-momentum tensor r^ v . The reason is as follows: The brane sources (in our model: r^ v and 
70 ) appear explicitly in the junction conditions only, which determine the jumps in the values of the 
Lie derivatives (in our model: the extrinsic curvature, K^, and the Lie derivative of the dilaton, 
£ n 4>) of the brane fields (the induced metric tensor and the dilaton field 0, respectively). The 
left hand sides of the junction conditions ( [241 and ( [251) are given by polynomials of the order 
(2A r max — 1) with only odd powers of the Lie derivatives. Although in general there can be several 
different solutions of these equations, we are interested only in those which vanish in the limit of 
vanishing sources (in the absence of sources the jumps in the values of the Lie derivatives have to 
vanish). Such solutions can be written as series in r^ u and with vanishing constant terms. On 
the other hand, only even powers of the Lie derivatives (K^ u and £ n <fr) are present in the (bulk) 
equations of motion. Consequently, in our effective brane equation ( [413]) the terms involving the 
sources are at least quadratic (or bilinear) in t^ v and r^. Moreover, among the bilinear terms there 
is no term proportional to t^t^ v . As we will show in the next section, such term is absent in the case 
of N m3uX = 1. It cannot appear for higher iV max as well, as the higher order corrections can change 
only those terms in the junction conditions solutions which are of higher order in the sources (we 
recall that one should consider only such solutions which have no constant terms when expanded 
in the sources). 

In the standard Einstein equation the energy-momentum tensor appears linearly only. In order 
to have such a term in our effective brane equations of motion (j3HJ) and P9|) we have to rewrite 
as a sum of the energy-momentum tensor t^ v associated with the fields we are interested in (e.g. the 
Standard Model fields which are usually assumed to be localized on a brane) and some "cosmological 
constant" A term: t^ w = t^ w + h^ u \. The result of such a redefinition will be discussed in more 
detail in the next section, which is devoted to the dilaton gravity with N max = 1 and iV max = 2. 



5 Examples 

The results presented in the previous sections are valid for models with corrections of arbitrarily 
high orders. In this section we write down those results explicitly^] for two simplest cases of 
-?Vm ax = 1 and iV max = 2. Although models with higher order corrections are the main topic of 
this work, we nevertheless want to present the results also for the lowest order theory. There are 
two reasons. First: the complexity of the calculations grows rapidly with the order of corrections. 



10 Due to their complexity, explicit formulae leading to the -/V max = 2 Einstein-like brane equation are moved to 
the appendix. 



19 



Thus, it is more instructive to discuss the main features of our procedure and results in the simplest 
situation with iV max = 1. Second: we have obtained new results even for the lowest order theory. 
In our approach we treat the dilaton field on the same footing as the metric tensor, which has not 
been done before. The results on the effective brane equations for the dilaton gravity obtained so 
far in the literature were not fully satisfactory^. 

In the present work we give only the general formulae for -/V max = 1 and iV max = 2. Presentation 
and discussion of some specific examples are postponed to a future publication. 

5 . 1 N max — 1 

Let us illustrate the main features of the construction of the effective brane equations of motion by 
considering a simple example with iV max = 1. In this case the "+" and the "— " sides limits of the 
bulk tensor equation of motion projected on the brane (1301) have the form of 

[R^ + (DD)^ - i V [ R + 2(DD)((> - (D0) 2 } + «rV(0) V} + { 2 (KK)^ - K^(K - £ n <f>) 

-\h^(?>(KK) -{K- £ n <j>) 2 - 2a A V A </>) } - [£ n K^ - h^ v {(h£ n K) - £ n 2 <f) } 

whereas the "+" and the "— " sides limits of the dilaton equation of motion ( )33l) read 

R + 2(DD)(f) - (D(f)) 2 - 2aJ- 1 (V(<j>) - V'{<j>)) } + {3 (KK) - (K - £ n (f)) 2 - 2 a x V x 4 



= 0, (50) 

(+) 



-2{(h£ n K)- £ n 2 <j)} 



= 0. (51) 

-(+) 



Both of these equations are genuine bulk equations of motion. They just determine the second Lie 
derivatives of the fields {£ n K^ v and £ n 2 (j> in the third curly bracket in each equation) in terms of the 
fields (the metric tensor and the dilaton field <p, together with their derivatives along the brane 
directions, in the first curly bracket) and their first Lie derivative^ (K^ and £ n (p in the second 
curly bracket). The values of all the quantities in the first two curly brackets in each equation are 



11 For example, in ref. |14| not all the calculations were carried out in a fully covariant way. Terms containing the 
combination <z a Va</> were removed by a gauge choice. Moreover, the influence of the bulk scalar field on the brane 
dynamics was taken into account only in some approximation and its bulk behavior was not eliminated from the 

brane gravitational equations. 

12 The last term in the second curly bracket is a mixture of the first Lie derivative of the dilaton with its derivatives 
along the brane directions. Hence, a part of that term should be moved to the first curly bracket. However, this 
subtlety does not change any further reasoning. Moreover, as was already mentioned, this slightly problematic 
term, a x \7\(j>, does not appear in our final results due to the appropriately designed derivation of the effective brane 
equations. 



20 



just the Cauchy boundary conditions for the corresponding <i-dimensional second order differential 
equations. Without other equations or any additional assumptions about the bulk solution, these 
boundary conditions can be arbitrary. Thus eqs. (1501) and (15 ip do not give any constraints on the 
brane fields and </> in terms of the brane sources r^ and r^. To get any effective brane equations 
of motion we have to analyze the junction conditions and the bulk tensor equation of motion with 
at least one index contracted with that of the vector n M normal to the brane. 

The junction conditions ( 124]) and ( 125]) are very simple in the case of N mSiX = 1. Solving them 
we can express the jumps in the values of the extrinsic curvature, K^, and the Lie derivative of 
the scalar field, £ n (j), in terms of the brane sources, r^ u and r^, as 

[K(u] ± = «r X { h ^ r <t> ~ T ^y) , ( 52 ) 
[£ n (P) ± = a^((d-2)r 4> -r). (53) 

Subsequently, we consider two vector eqs. ( 134]) - one for each side of the brane. The difference 
of those two equations together with the junction conditions (152]) and ( 153]) gives the consistency 
condition (I46p for the sources, namely 

D%, + Dy ( - V,) = . (54) 

As the second combination of the "— " and the "+" side eqs. (I34p we take the "+" one. Hence, we 
obtain the following condition for the "+" side quantities: 

D» [K^] + - [K, u ] + -D V [K] + + D v [£ n <j>] + = . (55) 

Finally, we have to take into account two scalar eqs. ( 135]) with iV max = 1. The brane curvature 
R and the dilaton <f), together with the first and the second (covariant with respect to the brane 
metric h^ v ) derivatives of the latter all are continuous at the brane. Thus their contributions cancel 
in the difference of the "— " and the "+" side eqs. ( 135]) . Such a difference reduces to the equality 
= \{K — £ n (j)) 2 ] ± — [(KK)} ± , which, after employing the junction conditions (152]) and (153]) . can 
be rewritten as 

r"" [K^] + - [£ n( f>] + + \a^ 1 ((rr) - 2rr + (d - 2)r}) = , (56) 

where (rr) = r^r^. The "+" side eq. (135]) reads 

R + 2(DD)cf> - (D<P) 2 - 2a^V(<f>) = ([K] + - \£ n <P} + ) 2 - \{KK)} + . (57) 

The last six eqs. (|5^]) - (l5T]) are the only equations which can yield the effective brane equations of 
motion in a general case. The junction conditions (15^]) and (133]) determine the jumps in the values of 
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the Lie derivatives of the metric tensor and the dilaton field: [i^Ji and [£ n 0]±, respectively. These 
quantities do not appear in any of the remaining four equations. Equation (154"]) is a consistency 
condition on the brane sources t^ u and - corresponding to the covariant conservation of the 
energy momentum tensor in the standard theory of gravity. The remaining equations establish 
conditions for the directional limits of the values of the first derivatives (normal to the brane) of 
the brane fields - evaluated at one of the "sides of the brane" , which we chose to be the "+" side. 
Specifically, eq. ( 1561) relates and [,£„</>]+ to the brane sources r^ v and r^, while eqs. ( 1551) 

and (1571) relate them to the brane fields h^ u and 0. Equations (l55l) -()571) provide us with (d + 1) 
relations, i.e. less than the number of the (gauge-independent) degrees of freedom in [iiT M i/]+ and 
[£ n (f)]+. Thus, in general, eqs. fl55|) -f )57|) can be solved for arbitrary brane fields and sources. The 
junction conditions ( 15l2"j) and (155j) do not change this situation, as without any assumptions on 
the bulk solution (as e.g. the already mentioned and usually employed bulk Z 2 symmetry) the 
quantities [-R^]±, [£ n <P}±, [K fJjU ] + and [£ n <f)] + are all independent. 

The only equation which involves the brane quantities exclusively is given by the formula (154)) . 
However, it is usually considered as a consistency condition on the brane sources and not as a 
dynamical equation of motion. As we pointed out in the previous section, it can yield a dynamical 
equation if the brane Lagrangian Cb is complicated enough, e.g. when it involves brane localized 
kinetic terms for the bulk fields. 

Let us now check the implications of some usually assumed features of the bulk solution. The by 
far most popular assumption of this kind is the already discussed bulk Z 2 symmetry with the fixed 
point at the brane position. With this symmetry it is enough to consider eqs. (134"|) and fl35|) on one 
"side of the brane" only. The corresponding equations on the other side of the brane are fulfilled 
automatically. Hence, in the case of iV max = 1 we are left with eqs. (p)2]) -( l5~4"|) and fl57|) . Employing 
[i^i/]± = 2[K^ V } + and [»£ n 0]± = 2[£ n (p} + relations, which are due to the bulk Z 2 symmetry, together 
with the junction conditions (152]) and (153]) . we can rewrite eq. (|57|) in the following form: 

R + 2{DD)(j) - (D<P) 2 - 2a^ l V{(j)) = \af [ -(rr) + 2 tt> - (d - 2) rj ] . (58) 

It determines the brane curvature scalar R in terms of the brane sources and (the dilaton 
field derivative can be obtained from the consistency condition (l5l| ). 

We shall now construct explicitly the brane tensor equation of motion for iV max = 1 - following 
the procedure for arbitrary iV max , described at the end of subsection 14.21 We start with eliminating 
the second derivatives normal to the brane, £ n K fJ-l/ and £ n 2 (f), from the bulk tensor equation f )50|) 
and the bulk dilaton equation of motion ( 15T|) . The parameters b m and B m , defined in (|4T|) and 
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(1421) . can be easily calculated: 

i-iy. 

-l-m) 



«i T^TT , (59) 



B m = ai {l(d-3) m [R-K 2 + (KK)] + (d-2) m [{DD)<j> + K£ T 

-\{d - l) m [(D(j)) 2 + (£ n 0) 2 ] } - (d - l) m V + (1 - m) V . (60) 



Employing the above explicit formulae on the parameters b m and B m , and the definitions f[3"Tj) 
and flUJJ), with and £ n given by their limits and [£ n 0] + obtained from the junction 

conditions (152"}) and fl53l) supplemented by the relations fH7|) due to the bulk Z 2 symmetry, the 
brane tensor equation (H9l) reduces to 



■K 2 



d-2 TT V ( rr )tiu + h^ u ( 2 (tt) (d-l)(d-2) T 



d-3 , (d-2)(d-3) 2 

' d-1 ' '<$> 2(d— 1) '0 



(61) 



This equations is truly of the form of the (d— l)-dimensional Einstein equation with some corrections 
(which was not apparent for the general iV max formula f H9|) ). Three specific types of contributions 
can be discerned on its r.h.s. There are terms with the explicit ^-dependence, typical of the gravity 
theories with scalar fields. Moreover, the tensor E^ v represents the bulk influence on the dynamics 
at the brane. Finally, the last square bracket contains the contributions from the brane sources 
t^v and T0. These contributions are quadratic in the brane energy-momentum tensor (which 
is typical of the brane models) and its dilaton counterpart r^. In order to have terms linear in 
some energy-momentum tensor (as is the case in the standard Einstein gravity) we have to rewrite 
t^v as the already mentioned sum of the energy-momentum tensor r Mi , associated with the fields 
we are interested in and some "cosmological constant" A term: r^ v = T yw + h^ v \. With such a 
decomposition we get the Einstein-like brane equation of motion as 

R,u - \h^R=MGr^ - ^[(DD)^ - h, u (DD)<f)} - ^h^D^f + a^V^)] - E, u 



+K 2 



where we introduced 



d-2 TT V ( TT )^ + 2^V" ( TT ) (d-l){d-2)^^ v T (d-2 ^ d-1 T <p) T 

A 2 - (d - 3) A r* + rf) hj , (62) 



G= Z y ;Z 2 i (63) 



jd-3) A 
32(d-2)iral 

which can be interpreted as the effective brane Newton's constant. The contributions to eq. (16"2|) 
which are proportional to the brane metric, h^, and depend neither on the scalar field <fr, nor on 
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the tensor t^, should be interpreted as the effective brane cosmo logical constant: 

A = fcf ar V|, =0 - i « r 2 (¥ A 2 - (d - 3) A r,|, =0 + r||, =0 ) . (64) 

It depends on the brane Lagrangian Cb (via r^) and on the way in which we divide the brane 
energy-momentum tensor r pv into its "standard" part T pv and the "cosmological constant" term 
hf^X. In addition, a part of E pv proportional to h pu can be also treated as a contribution to the 
effective brane cosmological constant. 

The effective brane equation of motion (162]) has the same tensor character as the standard 
Einstein equation. This certainly does not mean that we found another bulk-independent effective 
brane equation of motion - in addition to (1581) . It is just a convenient way to parameterize the bulk 
influence by a single geometric quantity: the projection of the bulk Weyl tensor on the brane, E pu . 
Observe that there is no additional bulk influence due to the presence of the dilaton fields. 



5.2 



iVmax = 2 



The formulae for arbitrary order of corrections given in sections [3] and H] can be employed to obtain 
more explicit equations for any given N max . However, the complexity of the resulting expressions 
grows rapidly with A max . The effective brane equations of motion become quite intricate already 
for the -/V" max = 2 case. Nevertheless, they can be still obtained even in the most general case, 
i.e. without any additional assumptions on the bulk background. We impose only the usual Z 2 
symmetry. As was already underlined, the A max = 2 case is equivalent to the appropriate subset 
of higher order interaction terms of the effective action derived from string theories. 
In the A max = 2 case the junction conditions ([24]) and 025p take the following form: 



a x [h pv (K - £ n (f>) - Kfj, v ] 

- 2 a 2 { [V(# " £n<f>) ~ [R-K 2 + (KK) + 2 (DD)<P + 2 K£ n cf> - (D<j)) 2 
- 2 V K pa [R pa + (KK) pa + (DD) pa <f\ —2(K— £ n< j>) [R^ + (DD)^<j> -\ 
+ 2 K, p [RP + {KKX + (DD) p (f\ + 2 K vp [RP + (KK); + [DD\ 
+ 2 K0° [R wvu - K w K va \ + | [h pu {K - £ n <pf + 2{KKK)} 



(K - £ n (j)) + 2 a 2 {(K - £ n (j)) [R - K 2 + (KK) + 2 {DD)(f> + 2 K£ n (f) - (D<j>f 
+ \{K- £ n< j>f - 2K pu [RP» + \{KKY U + (DD)^ 



- (£n<f>) 2 } 
{KK) pu \ 

(65) 

(£n^) 2 ] 

(66) 



This differs from the results previously presented in the literature. For example, it is claimed in |15j that in the 
effective brane equations of motion it is easier to remove the dependence on the projected bulk Weyl tensor than 
the dependence on the bulk dilaton field. Our analysis clearly indicates the opposite. 
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where (KKK) denotes the trace of the third power of the extrinsic curvature: K u p K^K p x . There 
are two new features as compared to the lowest order case. First: the junction conditions (165)) 
and (166 p are no longer linear in [iC^] and [£ n <f)} + . They are now third order equations for these 
quantities and their tensor structure is much more complicated. Thus, obtaining an explicit result 
is considerably more difficult. Second: solving these junction conditions (explicitly or not) yields 
the jumps in the values of K^ v and £ n (f> as functions not only of the brane sources r pv and r^, but 
also of the brane curvature R and the dilaton 0. 

The simplest form of the effective brane equation of motion is given by the scalar eq. ( 1351) . In 
the model with the bulk Z 2 symmetry and iV max = 2 it reads: 

V -\ ax [R - K 2 + (KK) + 2 (DD)<f> + 2 K£ n <\> - (D(f)) 2 - (£ n <j)) 2 ] 
-\ a 2 { [R - K 2 + (KK) + 2 (DD)<j> + 2 K £ n 4> - (D4>) 2 - (£ n 4>) 2 ] 2 
- 4 [R^ - KK^ + (KK)^ + (DD)^ + K^£ n cf\ ■ 
■ [R^ - KK^ + (KKY U + (DDY v (j) + K» v £ n (f\ 
+ [iW - K W K VCJ + K^K vp ] [RTf" - K w K ua + K^K V "] } = , (67) 

with K^v (and all its contractions) and £ n cf) replaced by their "next to the brane" values + 
and [£ n (p} + , given by the solutions of the junction conditions ( l65l) and ( |66l) . This is the only bulk- 
independent brane equation of motion in addition to the consistency condition ( BET) on the brane 
sources and r^. It is also possible to derive the Einstein-like effective equation of motion, as 
discussed at the end of subsection 14.21 However, the result is truly complicated and we postpone 
presenting the appropriate formulae to the appendix. Those formulae indicate that employing the 
brane Einstein-like effective equation for a general case is rather problematic. However, it can be 
significantly simplified if we restrict our considerations to the situations with sufficient symmetries. 
This is quite typical of all theories of gravity. Usually only solutions with some specific symmetry 
properties are looked for. An alternative way to investigate some highly symmetric solutions is to 
use the simplest form of the effective brane equation (1BT|) instead of the Einstein-like form discussed 
in the appendix. 

6 Conclusions 

The starting point of the present analysis has been given by the (i-dimensional higher order dilaton 
gravity constructed previously [7] as a generalization of the Einstein-Lovelock theory and remaining 
in the close relation to the effective action in string theories if restricted to the gravity and the 
dilaton field. It was supplemented by a co-dimension 1 brane with general brane localized inter- 
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actions Lb included into the Lagrangian ([T]). The effective brane equations of motion for such a 
theory were constructed and discussed. All calculations were performed in the covariant approach. 

In order to obtain the effective brane equations of motion one has to start from the full bulk 
equations of motion (TH|) . (fT5|) and attempt to eliminate all quantities evaluated away from the 
brane position. In general this is not possible and the effective brane equations do not form a 
closed system. The dynamics of the brane fields (the induced brane metric tensor) and <fi (the 
dilaton field) depends usually not only on the brane sources and defined in but also on 
the bulk gravity solution. 

On the basis of the full bulk equations of motion derived from the Lagrangian density defining 
our model, we can obtain three types of equations involving fields determined either on or infinites- 
imally close to the brane. These equations represent the junction conditions and two directional 
limits of the bulk equations established when the brane is approached from the "+" and the "— " 
sides, respectively. The junction conditions (1241) and (1251) relate the brane fields, and <fr, and the 
sources, r MV and r^, to the across-the-brane jumps in the values of the extrinsic curvature and the 
Lie derivative (along the vector field n M ortho normal to the brane) of the scalar field: [^K"^^] _j_ and 
[£ n 4>} ± , respectively. The directional limits of the bulk equations (given by eqs. (1501) . ( 155]) - (155]) ) 
involve the brane fields (h^ u and 0), together with their first [K^v and £ n (p) and second (£ n K PiU 
and £n<$) Lie derivatives, but do not include the brane sources r MJ , and r^. An important point 
should be underlined: the number of the d- dimensional bulk gravitational equations of motion is 
bigger than the number of equations necessary for a (d — l)-dimensional brane gravity. These extra 
equations (corresponding to the consistency conditions for a Cauchy problem with the boundary 
conditions defined infinitesimally close to the brane) play a crucial role, as they do not depend on 
the second Lie derivatives of the fields, i.e. £ n K ptV and £ n 2 (j). 

Genuine effective brane equations of motion should relate the values of the fields, h^ v and 0, 
evaluated on the brane, to the brane sources, and r^. In addition, they should not depend on 
the bulk configuration. The crucial question is: How many such brane equations can be obtained 
by combining information from all bulk equations mentioned in the previous paragraph? The 
answer depends strongly on the symmetries we assume for the bulk theory and for the bulk and/or 
brane solutions. The minimal number of the effective brane equations is obtained when no such 
symmetries are assumed. The only effective brane equation in such a case is given by formula (1461) . 
However, this equation is usually treated not as a dynamical brane equation of motion, but rather as 
a consistency condition on the brane sources ([7]). In the pure gravity case (i.e. without the dilaton) 
it reduces simply to the covariant conservation of the brane localized energy-momentum tensor. 
However, there are models in which this "consistency" condition yields a dynamical equation of 
motion. This is the case when e.g. brane localized kinetic terms are present in the brane localized 
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Lagrangian Lb- 

The number of the effective brane equations of motion increases when we restrict the model 
by imposing some symmetries. The bulk Z 2 symmetry with the fixed point coinciding with the 
brane position is particularly frequently employed. In such additional effective brane 

equation of motion fl48jl appears. It is obtained from the bulk tensor equation of motion if both 
indices are contracted with the vector field n M normal to the brane - after applying the appropriate 
junction conditions. This is the only effective brane equation of motion in models for which the 
"consistency" condition (1461) is not dynamical. The importance of this equation depends on the 
class of solutions we are interested in. For example, it is all we need when considering maximally 
symmetric brane solutions. 

In many models the number of the effective brane equations of motion is (much) smaller than 
the number of independent components of the Einstein equation in a (d — l)-dimensional space- 
time. Nevertheless, it is useful to derive an Einstein-like (tensor) brane equation. Such an equation 
for the higher (arbitrary) order dilaton gravity with the bulk Z 2 symmetry is given by eq. (149]) . 
The entire dependence of the brane dynamics on the bulk gravity solution is encoded in the Weyl 
tensor projected on the brane, E^ v , which appears in eq. (149 p implicitly via the parameters defined 
by formulae (142]) and (|43|) . It should be stressed that, contrary to some previous claims, those 
effective brane gravitational equations do not depend on the bulk scalar solution. 

Our general results obtained for corrections up to order 2A max (arbitrary, as long as it is not 
higher than the space-time dimensionality) in derivatives are presented in a very compact notation 
based on the generalizations (T5]) and ([8]) of the trace operation. Those results are rewritten explicitly 
in the conventional notation for the two simplest cases of N max = 1 and A max = 2. Although the 
analysis of the higher order theories is the main topic of our work, two reasons motivated us to 
address also the lowest order theory (i.e. with A max = 1), which has been already considered by 
other authors. First, it is useful as the simplest illustration of our general, non-trivial procedure. 
Second, we improved the analysis presented so far in the literature even for this lowest order theory. 
The resulting explicit formulae for the effective brane Newton's and cosmological constants are given 
by eqs. (IB5|) and (El]), respectively. 

The case of A max = 2 is the Einstein-Gauss-Bonnet gravity interacting with a scalar field 
(which is also self-interacting) via terms with up to four derivatives. Similarly to the A max = 1 
case, the effective brane equations of motion do not involve quantities dependent on the bulk 
scalar solution. The total bulk influence enters again through the projected Weyl tensor E^ u . The 
effective brane equations of motion for such a theory have not been presented before. They are quite 
lengthy when the generalized traces are explicitly calculated (all necessary formulae are collected 
in the appendix). Although the effective brane equations are rather complicated for a general case, 
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they simplify substantially if highly symmetric branes are considered. Applications of the derived 
equations for such symmetric models are postponed to a future publication. 
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Appendix 

In the appendix we collect the formulae appearing in the Einstein-like effective brane equation fl49|) 
for the case of iV max = 2. This equation, after writing down explicitly the sum over N, takes the 
form of 

(d-3)(2rf-3) 7 leiYf (M) + —h fl h' y T (V) + 4=$h T(M) 

(d-l)(d-2) U P aV 2 "n d-2 lb p lb a 1 P u \ JVL ) ^ d-3 P o ' P v \ ' I ^ d~l lb P a 1 \ JVl ) 

i of j o^ biBo-boBi 2 6igi - b 2 B 

+ 6 > b b 2 -bl V+ d - 2 b b 2 -bl V 

(M 2 ) + ~^KKT, U (VM) + gVT (M 2 ) - ^KKT, V (A/") 

+ ^ b Z~-T KKT ^ {KM) + A ^ b th~-t KKT ^ {M) } =0 - (A - l} 

Performing the summations over N in definitions ( 14 ip and (1421) . we obtain the following expressions 
for the parameters b m and B m : 

b = ai + 2a 2 T(M), (A.2) 

b 1 =a 1 (d-l) + 2a 2 T(KM), (A.3) 

b 2 =a 1 (d-2)(d-l) + 2a 2 T ((K) 2 M) , (A.4) 

B =-V + V + ai T(M) + a 2 [T(M 2 ) + ^T(VM) -8T(A")] , (A.5) 

Bl = -(d - i)v + ai r (km) + a 2 [r {km 2 ) + (kvm) - s r (KAf) ] . (a.6) 
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Various generalized traces present in the above formulae should be replaced with the following 
explicit expressions: 

T(M) = R- A 2 + (KK) + 2(DD)(j) + 2K£ n (f) - (Dcj)) 2 - (£ n 0) 2 , (A.7) 

T(M 2 ) = [R-K 2 + (KK) + 2{DD)(f) + 2K£ n (j)-(D(t)) 2 -(£ n <t)) 2 ] 2 

- 4 [R^ - KK^ + (KK)fj, v + (DD)^<f> + K^£ n 

■ [R^ - KK» V + [KKY V + (OOf > + K^£ n 
+ [R pupa - K W K V<T + K^K vp ] [R^ - K w K va + K^K V "] , (A.8) 

T (VM) = -± l [R-K 2 + {KK)} [R - K 2 + (KK) + {DD)<f> + K£ n <f] 

- 2 [R pv - KK pu + (KK) pu + (DD) pu (p + K pv £ n (f] 

■ [R^ - KK» V + (KKY V + (d - 2)E" V ] , (A.9) 

T (AT) = [D p K - D p K? + K£D p <j> - D^£ n <j>] [D»K - D°K» + K^<j> - D»£ n <f] 

- (D p K° p )(D»K?) + {D,K;){D?KZ) , (A. 10) 

T M) = (d - 3) [i? - K 2 + (ATA)] + 2{d - 2) [{DD)<f> + K£ n (f] 

-(d-l)[(D(f>) 2 + (£ n( f>) 2 ], (A.ll) 

T (/£ M 2 ) = (d - 5) [R - K 2 + (M)] 2 + 2(d - 4) [R - K 2 + (M)] [(,0,0)0 + A£ n 0] 
+ (d - 3) (4 [(DD)4> + K£ n <p] 2 -[R-K 2 + (KK)] [(O0) 2 + (£ n <P) 2 ] ) 

- 4(d - 2) [(OO)0 + A£ n 0] [(D0) 2 + (£ n 0) 2 ] + (d - 1) [(D<P) 2 + (£ n 0) 2 ] 2 

- 4(d - 5) [i^ - KK pv + [i^ - KK^ U + (M)"^] 

- 8(d - 4) [i^ - KK pv + (AA%,J [(^^) M > + A0] 

- 4(d - 3) [(DD)^ + A M „£ n 0] [(55)^ + K^£ n( f\ 

+ (d - 5) [i^ - K pp K va + A^ A^J - K W K V(T + A^ A^>] , (A. 12) 

T (/it PM) = [i? - A 2 + (KK)] 2 + [R - K 2 + (AA)] [(£>£>)<£ + A£ 

- 2 ((d - 4) [R^ - KK pv + (AA)^] + (d - 3) [(DD)^ + K jlv £ n 

■ \BT - KK^ + (KKy v + (d- 2)E^] , (A. 13) 

T (KM) = (d-4)([D p K-D p K?] [D*K-D°K£]-(D^)(D*Ka + (D,jq)(D>KSj) 
+ 2(d-3)K;[(D u K)(D^) - (D P K^)(D^)] 

+ (d - 2) [KPD p <f) - D p £ n <f)] [K^D'cf) - D»£ n <j>] , (A.14) 
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T((K) 2 M) = (d-4)(d-3)[R-K 2 + (KK)]+2(d-3)(d-2)[(DD)(j) + K£ n 

-(d-2)(d- 1) [Wf + (£ n 0) 2 ] , (A.15) 

h^KT^iM) = h pa T(M)-2[R pa -KK pa + (KK) p(T + (DD) pa( f> + K p<T £ n( f>], (A.16) 

K h «T^ (V) = [R-K 2 + (KK)] - [R pa - KK pa + (KK) pa ] , (A. 17) 

W p K 7% {M 2 ) = h pa T {M 2 ) - 4 [R pa - KK pa + {KK) pa + (DD) p(J + K pa £ n (f)] T (M) 
+ 8 [R pp - KK pp + (AA) W + (DD) pp( f> + J^£„0] 

■ [R\ - KK\ + {KKy G + + tf'*^] 

+ 8 \Rpanv - K pp K av + if^A.J \BT - KK^ + (if Kf + + K^£ n cf] 

- 4 [i^ A - K pv K p \ + if pA if M J [i^ - K a v K " A + if CT A A , (A. 18) 

h»KT pu (N) = h p(T T(N) + 2{D p K v p )(D> 1 K U(T ) 

- 2 [D p K - D V K;] [D»K pa - D p K»] + [D p K; - D p K;] [D a K» - D U K£] 

- [D p K - D p K» + K^D p <f) - D p £ n cf] [D a K - D U K» + K v „D v< t> - D a £ n <j>] , (A.19) 

W p h v a T pv (KM) = V (i d ~ 4 ) [ R ~ K 2 + (KK)] + 2(d - 3) [(£>£>)<£ + A£ n 0] 

-(rf-2)[( J D0) 2 + (i? n ,0) 2 ]) 

- 2 (d - 3) [i? pCT - AA pff + (KK) pa ] - 2(d - 2) [(DD)^ + K pa £ n <P] , (A.20) 

frjX 7% (VM) = h pa T (VM) - ^ p K T pv (KM) - [R pa - KK P<J + (KK) pa + (d - 2)E pa ] 
-2[R-K 2 + (KK)] [R pa - KK pa + (KK) pa + (DD) pa <f) + K pa £ n cf] 
+ 2 [R PI1 - KK pp + (if K) PI1 + (d- 2)E pp ] [R£ - KK» + {KK)£ + (DD)£</> + K^£ n 
+ 2[R£- KK% + (KKY a + (d- 2)E£\ [R pp - KK PP + (KK) PP + (DD) PP + A PM £ n 
+ 2 [i^ - if^A^ + K pv K ap \ [FT - KK^ + (KK)^ + (d - 2)E flu ] . (A.21) 

By substituting eqs. ( 1A.2[) -( |A~2~T!) into ( lA.lj) and employing the solutions of the junction conditions 
(1651) and (1661) . the effective brane equation of motion in the Einstein-like form can be obtained. 
However, the result is quite intricate and will not be given here. Moreover, the formulae collected 
in this appendix should be compared with the relatively simple eq. f )67|) . This shows how big price, 
in terms of complication, has to be paid in order to derive the effective brane equation of motion 
in the Einstein-like form, instead of confining to eq. (]67p . 

Similarly to the A max = 1 case, the whole bulk dependence of the brane dynamics is encoded 
in the projected Weyl tensor E pv . However, this dependence is quite complicated for A max = 2. 
Specifically, E pv enters the brane Einstein-like equation of motion OA.ip through the generalized 
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trace h)th v a T^u (VM.), appearing in eq. (IA.1|) and given by (IA.21j) . as well as through the generalized 



traces T(VM) and T{h*JPM), present in the definitions of B Q and B\ and given by (1A.9|) and 
(IA. 13j) . respectively. Although the projected Weyl tensor E^ u enters these formulae only linearly, 
it is involved in intricate contractions with other tensors. Furthermore, these tensors contain the 
extrinsic curvature K^ u , i.e. a complicated solution of the junction conditions f )65|) and fl66|) - which 
in general cannot be solved explicitly. 

As for N max > 2, it is in principle possible to write down the Einstein-like effective brane 
equation of motion explicitly, but they become practically intractable. 
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